Cold collisions of OH and Rb. I: the free collision 



Manuel Lara/'E John L. Bohn, 2 >0 Daniel E. Potter, 3 Pavel Soldan, 4 '[3 and Jeremy M. Hutson 3 'H 

'JILA, University of Colorado, Boulder, CO 80309-0440 
2 JILA and Department of Physics, University of Colorado, Boulder, CO 80309-0440 
3 Department of Chemistry, University of Durham, England 
Department of Chemistry, Faculty of Nuclear Science and Physical Engineering, Czech Technical University, Czech Republic 

(Dated: February 2, 2008) 

We have calculated elastic and state-resolved inelastic cross sections for cold and ultracold colli- 
sions in the Rb( 1 S) + OH( 2 n 3 / 2 ) system, including fine-structure and hyperfine effects. We have 
developed a new set of five potential energy surfaces for Rb-OH( 2 II) from high-level ab initio elec- 
tronic structure calculations, which exhibit conical intersections between covalent and ion-pair states. 
The surfaces are transformed to a quasidiabatic representation. The collision problem is expanded in 
a set of channels suitable for handling the system in the presence of electric and/or magnetic fields, 
although we consider the zero-field limit in this work. Because of the large number of scattering 
channels involved, we propose and make use of suitable approximations. To account for the hyper- 
fine structure of both collision partners in the short-range region we develop a frame-transformation 
procedure which includes most of the hyperfine Hamiltonian. Scattering cross sections on the order 
of fO -13 cm 2 are predicted for temperatures typical of Stark decelerators. We also conclude that 
spin orientation of the partners is completely disrupted during the collision. Implications for both 
sympathetic cooling of OH molecules in an environment of ultracold Rb atoms and experimental 
observability of the collisions are discussed. 

PACS numbers: 34.20.Mq,34.50.-s,33.80.Ps 



I. INTRODUCTION 

The possibility of producing translationally ultracold 
molecules has recently generated great anticipation in the 
field of molecular dynamics. Attractive applications in- 
clude the possibility of testing fundamental symmetries 
, the potential of new phases of matter [3- C3- C3- 13| and 
the renewed quest for the control of chemical reactions 
through ultracold chemistry [E 0, HI • These endeavors 
are beginning to bear scientific fruit. For example, high- 
resolution spectroscopic measurements of translationally 
cold samples of OH should allow improved astrophysical 
tests of the variation with time of the fine-structure con- 
stant . The recent experimental advances 0, 0, 
have made theoretical studies of the collisional behavior 
of cold molecules essential [lj, L!-4 > both to interpret 
the data and to suggest future directions. 

Several approaches have produced cold neutral 
molecules to date, many of which are described in Ref. 
[Til The methods available can be classified into di- 
rect methods, based on cooling of preexisting molecules, 
and indirect methods, which create ultracold molecules 
from ultracold atoms. Among the direct methods, Stark 
deceleration of dipolar molecules in a supersonic beam 
[lii IttL ITsf and helium buffer-gas cooling are cur- 
rently leading the way. They reach temperatures of the 
order of 10 mK, and there are a wide variety of pro- 
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posals on how to bridge the temperature gap to below 
1 mK. These include evaporative cooling and even direct 
laser cooling. The idea of sympathetic cooling, where a 
hot species is cooled via collisions with a cold one, also 
seems very attractive and is being pursued by several 
experimental groups. Sympathetic cooling is a form of 
collisional cooling, which works for multiple degrees of 
freedom simultaneously, ft does not rely on specific tran- 
sitions, which makes if suitable for cooling molecules. 
Collisional cooling is also the basis for helium buffer-gas 
cooling. 

Sympathetic cooling of trapped ions has already 
been demonstrated |2fJ, using a different laser-cooled 
ionic species as the refrigerant. Cooling of polyatomic 
molecules to sub-kelvin temperatures with ions has also 
been reported. This technique is expected to be suit- 
able for cooling molecules of very high mass, including 
those of biological relevance [2l| . But the ease with which 
alkali metal atoms can be cooled to ultracold tempera- 
tures makes them good candidates to use as a thermal 
reservoir to cool other species. They have already been 
used to cool "more difficult" atomic alkali metal partners. 
For example, BEC for 41 K was achieved by sympathetic 
cooling of potassium atoms with Rb atoms |22|. A the- 
oretical study of the viability of this cooling technique 
for molecules is desirable. There have been a number 
of theoretical studies of collisions of molecules with He 
El 0, |H HI H| , in support of buffer-gas cooling, 
but only a few with alkali metals pjj, l30l l3ll l32t l3o| . To 
our knowledge, no such study has included the effects of 
hyperfine structure. 

The main objective of the present work is to study cold 
collisions of OH with trapped Rb atoms. OH has been 



2 



successfully slowed by Stark deceleration in at least two 
laboratories [!3,|3J] To cool the molecules further, sym- 
pathetic cooling by thermal contact with 87 Rb is an at- 
tractive possibility. Rb is easily cooled and trapped in co- 
pious quantities and can be considered the workhorse for 
experiments on cold atoms. Temperatures below 100 fiK 
are reached even in a MOT environment (70 /iK using 
normal laser cooling and 7 (iK using techniques such as 
polarization gradient cooling). 

The cooling and lifetime of species in the trap depends 
largely on the ratio of elastic collision rates (which lead 
to thermalization of the sample) to inelastic ones. The 
latter can transfer molecules into non-trappable states 
and/or release kinetic energy, with resultant heating and 
trap loss. The characterization of the rates of both kinds 
of process is thus required. Since applied electric and 
magnetic fields offer the possibility of controlling colli- 
sions, it is very important to know the effects of such 
fields on the rates. At present, nothing is known about 
the low-temperature collision cross sections of Rb-OH or 
any similar system. 

Rb-OH can be considered as a benchmark system for 
the study of the feasibility of sympathetic cooling for 
molecules. Many molecule- alkali metal atom systems 
have deeply bound electronic states with ion-pair char- 
acter [3(| [331 and have collision rates that are enhanced 
by a "harpooning" mechanism. Both the atom and the 
diatom are open-shell doublet species, and can interact 
on two triplet and three singlet potential energy surfaces 
(PES). In addition, the OH radical has fine structure, 
including lambda-doubling, and both species possess nu- 
clear spins and hence hyperfine structure. Thus Rb-OH is 
considerably more complicated than other collision sys- 
tems that have been studied at low temperatures. In 
previous work we advanced the first estimates of cross 
sections (for both inelastic and elastic collisions), based 
on fully ab initio surfaces, for the collision of OH radicals 
with Rb atoms in the absence of external fields |33j • Here 
we provide details of the methodology used and discuss 
the potential surfaces and the state-resolved partial cross 
sections. 

This paper is organized as follows: Section II describes 
the calculation of ab initio PES for Rb-OH. Details of the 
electronic structure calculations are given and the meth- 
ods used for diabatization, interpolation and fitting are 
described. The general features of the resulting surfaces 
are analyzed. Section III describes the exact and approx- 
imate theoretical methodologies used for the dynamical 
calculations. Section IV presents the resulting cross sec- 
tions and discusses the possibility of sympathetic cooling. 
We also comment on the role expected for the harpoon- 
ing mechanism. Section 5 summarizes our results and de- 
scribes prospects for future work. Further details about 
electronic structure calculations and other channels basis 
sets to describe the dynamics are given in Appendixes 1 
and 2 respectively. 



II. POTENTIAL ENERGY SURFACES 

We have used ab initio electronic structure calculations 
to obtain PES for interaction of OH( 2 n) with Rb. The 
ground A 2 n state of OH has a 7r 3 configuration, while 
the ground 2 S state of Rb has 5s 1 . At long range, linear 
RbOH thus has *n and 3 n states. At nonlinear config- 
urations, 1 H splits into 1 A' and 1 A", with even and odd 
reflection symmetry in the molecular plane, whereas 3 n 
splits into 3 A 1 and 3 A". 

At shorter range, the situation is more complicated. 
The ion-pair threshold Rb + + OH~ lies only 2.35 eV 
above the neutral threshold. The corresponding 1 E + 
( A') ion-pair state drops very fast in energy with de- 
creasing distance because of the Coulomb attraction. Be- 
low, Jacobi coordinates (R, 9) will be used, R being the 
radial distance between the atom and the OH center of 
mass and 9 the angle between this line and the inter- 
nuclear axis. At the linear Rb-OH geometry, the ion- 
pair state crosses the covalent (non-ion-pair) state near 
R — 6.0 A, as shown in Figure ^ At nonlinear geome- 
tries, the ion-pair state has the same symmetry ( l A') as 
one of the covalent states, so there is an avoided cross- 
ing. There is thus a conical intersection between the two 
1 A' states at linear geometries, which may have major 
consequences for the scattering dynamics. 

The A' electronic wavefunctions near the conical in- 
tersection are made up from two quite different config- 
urations, so that a multiconfiguration electronic struc- 
ture approach is essential to describe them. We have 
therefore chosen to use MCSCF (multiconfiguration self- 
consistent field) calculations followed by MRCI (multircf- 
erence configuration interaction) calculations to charac- 
terize the surfaces. The electronic structure calculations 
initially produce adiabatic (Born-Oppenheimer) surfaces, 
but these are unsuitable for dynamics calculations both 
because they are difficult to interpolate (with derivative 
discontinuities at the conical intersections) and because 
there are nonadiabatic couplings between them that be- 
come infinite at the conical intersections. We have there- 
fore transformed the two 1 A' adiabatic surfaces that cross 
into a diabatic representation, where there are non-zero 
couplings between different surfaces but both the poten- 
tials and the couplings are smooth functions of coordi- 
nates. 

The electronic structure calculations are carried out 
using the MOLPRO package [38| . It was necessary to 
carry out an RHF (restricted Hartree-Fock) calculation 
to provide initial orbital guesses before an MCSCF calcu- 
lation. It is important that the Hartree-Fock calculation 
gives good orbitals for both the OH ir and Rb 5s orbitals 
at all geometries (even those inside the crossing, where 
the Rb 5s orbital is unoccupied in the ground state). 
In addition, it is important that the OH n orbitals are 
doubly occupied in the RHF calculations, as otherwise 
they are non-degenerate at linear geometries at the RHF 
level, and the MCSCF calculation is unable to recover 
the degeneracy. To ensure this, we begin with an RHF 
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FIG. 1: RbOH adiabatic potential curves from MRCI calcula- 
tions, showing crossing for the Rb-OH linear geometry (a) and 
avoided crossing for a slightly nonlinear geometry (9 = 159°, 
(b)). Panel (c) shows an expanded view of the curves in the 
center panel. 



calculation on RbOH rather than neutral RbOH. 

For Rb, we use the small-core quasirelativistic effective 
core potential (ECP) ECP28MWB [H| with the valence 
basis set from Ref. j40j. This treats the 4s, 4p and 5s elec- 
trons explicitly, but uses a pseudopotential to represent 
the core orbitals. For O and H, we use the aug-cc-pVTZ 
correlation-consistent basis sets of Dunning |4l| in un- 
contracted form. Electronic structure calculations were 
carried out at 275 geometries, constructed from all com- 
binations of 25 intermolecular distances R and 11 angles 9 
in Jacobi coordinates. The 25 distances were from 2.0 to 
6.0 A in steps of 0.25 A, from 6.0 to 9.0 A in steps of 0.5 A 
and from 9 A to 12 A in steps of 1 A. The OH bond length 
was fixed at r = 0.9706 A. The 11 angles were chosen to 
be Gauss-Lobatto quadrature points (4^] , which give op- 
timum quadratures to project out the Legendre compo- 
nents of the potential while retaining points at the two 
linear geometries. The linear points are essential to en- 
sure that the A' and A" surfaces are properly degenerate 
at linear geometries: if we used Gauss-Legendre points 
instead, the values of the A' and A" potentials at linear 
geometries would depend on extrapolation from nonlin- 
ear points and would be non-degenerate. The Gauss- 
Lobatto points correspond approximately to 6 = 0, 20.9, 
38.3, 55.6, 72.8, 90, 107.2, 124.4, 141.7, 159.1 and 180°, 
where 8 — is to the linear Rb-HO geometry. The cal- 
culations were in general carried out as angular scans at 
each distance, since this avoided most convergence prob- 
lems due to sudden changes in orbitals between geome- 
tries. 



A. Singlet states 

We carried out a state-averaged MCSCF calculation of 
the lowest 3 singlet states of neutral RbOH (two 1 A' and 
one 1 A"). Molecular orbital basis sets will be described 
using the notation (ua 1 , tia"), where the two integers in- 
dicate the number of A' and A" orbitals included. The 
orbital energies are shown schematically in Figure^ The 
MCSCF basis set includes a complete active space (CAS) 
constructed from the lowest (10,3) molecular orbitals, 
with the lowest (5,1) orbitals closed (doubly occupied in 
all configurations). The MCSCF calculation generates a 
common set of orbitals for the 3 states. The calculations 
were carried out in C s symmetry, but at linear geome- 
tries the two components of the n states are degenerate 
to within our precision (lO -8 ^). 

For cold molecule collisions, it is very important to 
have a good representation of the long-range forces. 
These include a large contribution from dispersion (in- 
termolecular correlation), so require a correlated treat- 
ment. We therefore use the MCSCF orbitals in an MRCI 
calculation, again of the lowest three electronic states. 
The MOLPRO package implements the "internally con- 
tracted" MRCI algorithm of Werner and Knowles [43| . 
The reference space in the MRCI is the same as the active 
space for the MCSCF, and single and double excitations 
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FIG. 2: Schematic representation of the RbOH molecular or- 
bitals from MCSCF calculations. 



are included from all orbitals except oxygen Is. As de- 
scribed in Appendix 1, the two 1 A 1 states are calculated 
in a single MRCI block, so that they share a common 
basis set. 

We encountered difficulties with non-degeneracy be- 
tween the two components of the 1 II states at linear ge- 
ometries. These are described in Appendix 1. However, 
using the basis sets and procedures described here, the 
non-degeneracies were never greater than 90 [iEh in the 
total energies for distances R > 2.25 A (and consider- 
ably less in the interaction energies around the linear 
minimum) . 



B. Transforming to a diabatic representation 

As described above, the two surfaces of 1 A' symmetry 
cross at conical intersections at linear geometries. For 
dynamical calculations, it is highly desirable to trans- 
form the adiabatic states into diabatic states (or, strictly, 



quasidiabatic states). MOLPRO contains code to carry 
out diabatization by maximizing overlap between the dia- 
batic states and those at a reference geometry. However, 
this did not work for our application because we were 
unable to find reference states that had enough overlap 
with the lowest adiabats at all geometries. We therefore 
adopted a different approach, based on matrix elements 
of angular momentum operators. We use a Cartesian co- 
ordinate system with the z axis along the OH bond. At 
any linear geometry, the n component of A 1 symmetry 
is uncontaminated by the ion-pair state, and the matrix 
elements of L z are 



CA'mL^A") = i; 
fA'^lL^A") = 0. 



(1) 



At nonlinear geometries, the actual X A' states can be 
represented approximately as a mixture of X and n com- 
ponents, 



cos <p sin < 
- sin d> cos < 



(2) 



where the "singlet" superscripts have been dropped to 
simplify notation. If Eq. [2] were exact, the matrix ele- 
ments of L z would be 



(1A'\L Z \A") = icos0; 
(2A'\L Z \A") = isinc/). 

The mixing angle <j) would thus be given by 



= tan 



■ i (2A'\L Z \A" 
(1A'\L Z \A" 



(3) 



(4) 



In the present work, we have taken the mixing angle to be 
defined by Eq.^J using matrix elements of L z calculated 
between the MRCI wavefunctions. This gives a mixing 
angle that, for linear geometries, is 4> — at long range 
and <f) = 7r/2 at short range (inside the crossing). 

One complication that arises here is that the signs of 
the three wavefunctions are arbitrary, and may change 
discontinuously from one geometry to another. The signs 
of the matrix elements obtained numerically by MOL- 
PRO are thus completely arbitrary. It was therefore nec- 
essary to pick a sign convention for the matrix elements 
at linear geometries and adjust the signs at other geome- 
tries to give a smoothly varying mixing angle. 

It should be noted that this diabatization procedure is 
not general, and will fail if there is any geometry where 
both the numerator and the denominator of Eq. 01 arc 
small. Fortunately, this was not encountered for RbOH. 
The sum of squares of the two matrix elements of L z 
was never less than 0.99 at distances from R — 3.0 A 
outwards, and never less than 0.7 even at R = 2.0 A. 

The mixing angles obtained for the singlet states of 
RbOH are shown as a contour plot in Figure [21 As ex- 
pected, 4> changes very suddenly from to 90° at lin- 
ear and near-linear geometries, but smoothly at strongly 
bent geometries. 
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FIG. 3: Contour plot of the diabatic mixing angle (j> (in de- 
grees ) for the 1 A' states of RbOH 

Once a smooth mixing angle has been determined, 
the diabatic potentials and coupling surface are obtained 
from 
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The diabatization was carried out using total elec- 
tronic energies (not interaction energies). The two sin- 
glet diabatic PES that correlate with Rb( 2 5) + OH( 2 n) 
are shown in Figure 01 and the diabatic ion-pair surface 
and the coupling potential are shown in Figure |SJ Some 
salient characteristics of the surfaces are given in Tabled 
As required, the 1 A" and 1 A' covalent states are degen- 
erate at the two linear geometries, with a relatively deep 
well (337 cm" 1 ) at Rb-OH and a much shallower well at 
Rb-HO. At bent geometries, it is notable that the poten- 
tial well is broader and considerably deeper for the 1 A" 
state than for the 1 A' state; indeed, the 1 A I state has a 
linear minimum, while the 1 A" has a bent minimum at 
9 = 128° with a well depth of 405 cm -1 . In spectroscopic 
terms, this corresponds to a Renner- Teller effect of type 
1(b ) in the classification of Pople and Longuet-Higgins 

The fact that the 1 A" state is deeper than the 1 A' state 
is somewhat unexpected, and will be discussed in more 
detail in the context of the triplet surfaces below. The 
fact that the surface for the covalent 1 A' state is slightly 
repulsive between R = 4 and 9 A and between 6 = 40° 
and 150° is an artifact of the diabatization procedure 
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FIG. 4: Contour plot of the diabatic 1 A' (panel (a)) and 1 A" 
(panel (b)) covalent (non- ion-pair) potential energy surfaces 
for RbOH from MRCI calculations. Contours are labeled in 
cm -1 . 



TABLE I: Characteristics of the diabatic PES for RbOH. 



^A" T A' T A n T ~A' T A' 

covalent ion-pair 
well depth/cm -1 615 511 405 337 26260 

distance of minimum/A 3.185 3.170 3.226 3.230 2.407 
angle of minimum/ ° 123 180 128 180 180 
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FIG. 5: Contour plot of the 1 A 1 diabatic ion-pair potential 
energy surface (panel (a)) and the diabatic coupling potential 
(panel (b)) for RbOH from MRCI calculations. Contours are 
labeled in cm - . 



and should not be given physical significance: the choice 
of mixing angle l@J is an approximation, and a slightly 
different choice would give different diabats and coupling- 
terms, but of course corresponding to the same adiabatic 
surfaces. 

The ion-pair state has a deep well at the RbOH geom- 
etry and a rather shallower one at RbHO, as expected 
from electrostatic considerations. The region around the 



minimum of this surface has been characterized in more 
detail by Lee and Wright p5|. 

It is notable that the coupling potential is quite large, 
peaking at 2578 cm" 1 at R = 3.87 A and 9 = 83°, and 
is thus larger than the interaction energy for the two 
covalent surfaces at most geometries. It may be seen from 
Figure |3| that the E — II mixing it induces is significant 
at most distances less than 8 A. 



C. Triplet states 

The triplet states of RbOH are considerably simpler 
than the singlet states, because there is no low-lying 
triplet ion-pair state. There is thus no conical intersec- 
tion, and diabatization is not needed. 

Single-reference calculations would in principle be ade- 
quate for the triplet states. Nevertheless, for consistency 
with the singlet surfaces, we carried out MCSCF and 
MRCI for the triplet surfaces as well. The (10,3) ac- 
tive and reference spaces were specified as for the singlet 
states, except that there is only one relevant 3 A 1 state and 
the state average in the MCSCF calculation is therefore 
over the lowest two states. 

Contour plots of the 3 A 1 and 3 A" surfaces are shown in 
Figure EJ As for the singlet states, it is notable that the 
3 A" state lies below the 3 A' state at nonlinear geometries. 
This ordering is different from that found for systems 
such as Ar-OH [H H3 and He-OH H^. In each case, 
the A' state corresponds to an atom approaching OH 
in the plane of the unpaired electron, while the A" state 
corresponds to an atom approaching OH out of the plane. 
For He-OH and Ar-OH, the 2 A' state is deeper than the 
2 A" state simply because there is slightly less repulsion 
due to a half-filled it orbital than due to a doubly filled 
7r orbital. Since these systems are dispersion-bound, and 
the dispersion coefficients are similar for the 2 A' and 2 A" 
states, the slightly reduced repulsion for the 2 A' state 
produces a larger well depth. 

Rb-OH is quite different. The long-range coefficients 
still provide a large part of the binding energy of the 
covalent states, but the equilibrium distances (around 3.2 
A, Tabled) are about 1 A shorter than would be expected 
from the sum of the Van der Waals radii of Rb (2.44 A) 
and OH (1.78 A, obtained from the Ar value of 1.88 A 
and the Ar-OH equilibrium distance of 3.67 A |47|V 

The qualitative explanation is that at nonlinear geome- 
tries there is significant overlap between the Rb 5s orbital 
and the OH 7r orbital of a' symmetry, forming weakly 
bonding and antibonding molecular orbitals (MOs) for 
the RbOH supermolecule. The effect of this is shown in 
Figured I n the A' states, the bonding and antibonding 
MOs are equally populated and there is no overall stabi- 
lization. However, for the A" states the bonding MO is 
doubly occupied and the antibonding MO is singly occu- 
pied. This gives a significant reduction in the repulsion 
compared to a simple overlap-based model. 

Even at linear Rb-OH geometries, the repulsion is re- 
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FIG. 6: Contour plot of the 3 A' (panel (a)) and 3 A" (panel 
(b)) PES for RbOH from MRCI calculations. Contours are 
labeled in cm - . 



duced by a similar effect involving the Rb 5s orbital and 
the highest occupied OH a orbital (6a' in Figure^. 
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FIG. 7: Orbital occupations for the covalent (non-ion-pair) 
states of RbOH at nonlinear geometries, 90° < 8 < 180°, 
showing the origin of the reduced repulsion for A" states com- 
pared to A' states. 



ing the counterpoise approach |49j • In addition, there are 
Rb 5p orbitals that lie 1.87 eV above the ground state. 
This is below the ion-pair energy at large R, so that a 
(10,3) active space includes different orbitals asymptot- 
ically and at short range. It is thus hard to calculate 
asymptotic energies that are consistent with the short- 
range energies directly. Nevertheless, for collision calcu- 
lations we ultimately need interaction energies, relative 
to the energy of free Rb( 2 S) + OH( 2 n). 

To circumvent this problem, we obtained angle- 
dependent long-range coefficients C§{9) and C 7 {9) for the 
triplet states of RbOH and used these to extrapolate from 
12 A outwards at each angle. We carried out RCCSD cal- 
culations (restricted coupled-cluster with single and dou- 
ble excitations) on the 3 A' and 3 A" states at all angles 
at distances of 15, 25 and 100 A. RCCSD was chosen in 
preference to RCCSD(T) for greater consistency with the 
MRCI calculations. Coupled cluster calculations are size- 
extensive, so in this case the interaction energies V(R, 9) 
for the 3 A 1 and 3 A" states were calculated including the 
counterpoise correction 49]. The interaction energy at 
100 A was found to be non-zero (about — 1.38 fJ.Eh), but 
was the same to within 10~ 9 Eh for both 3 A' and 3 A" 
states and at all angles. For each pair of angles 9 and 
7r — 9, the energies were fitted to the form 

E(R,0) = E^-C^R-* -C 7 (9)R- 7 -C 8 {9)R- S , (6) 



with the constraints that 



C 6 (tt-9) = C 6 (9) 
C 7 (ir-9) = -C 7 {9). 



(7) 

(8) 



Sums and differences of the long-range coefficients for the 
two states, 



D. Converting MRCI total energies to interaction 
energies 

The MRCI procedure is not size-extensive, so cannot 
be corrected for basis-set superposition error (BSSE) us- 



C° n (9) 
C 2 n (9) 



l r 

2 
1 



(9) 
(10) 



were then smoothed by fitting to the theoretical func- 
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TABLE II: Long-range coefficients for Rb-OH obtained by 
fitting to RCCSD calculations on triplet states. 



A 
Ct'/E h a b 325.0 
Cg 2 /Eha® 
C 7 °/Ef l a 
Cj 2 / Eh.a'a 



1035.4 



151.0 
1.9 

630.0 
-40.1 



results. Our approach to two-dimensional interpolation 
follows that of Meuwly and Hutson [5(j and Soldan et 
al. [SlJ. The interpolation was carried out first in R (for 
each surface and angular point) and then in 9. The in- 
terpolation in R used the RP-RKHS (reciprocal power 
reproducing kernel Hilbert space) procedure [35| with pa- 
rameters m = 5 and n — 3. This gives a potential with 
long-range form 

V q (R, 9) = -Cl(9)R- 6 - C q 7 {9)R- 7 - Cl(9)R- s . (17) 



tional forms 



C 6 » 



Ct{9) 



C°° + C 2 6 °P° (cos 9) 

cf Pi (cose) 

Cf P%{cos6) + Cf P 3 °(. 
CfP 3 2 (cos0), 



[cost 



(11) 

(12) 

(13) 
(14) 



where P£(cos9) are associated Legendre functions. The 
coefficients obtained by this procedure are summarized in 
Table|n| The resulting smoothed values of C%(d), C£(9), 
C°(9) and Cf(#) were then used to generate (9), 
Cf{9), C A '{9) and Cf(9). Finally, the MRCI total 
energies at R = 10 and 12 A for both singlet and triplet 
states were refitted to Eq. El with Cq(9) and CV(#) held 
constant at the smoothed values, to obtain an MRCI 
value of Eoo for each angle (and surface). These angle- 
dependent values of Poo were used to convert the MRCI 
total energies into interaction energies. 

It should be noted that the long-range coefficients in 
Rb-OH have substantial contributions from induction as 
well as dispersion. A simple dipole-induced dipole model 
gives C°° ind = C 2 ° nd = 133 E h a% and accounts for 40% 
of C°° and 90% of C 6 20 



E. Interpolation and fitting 



The values of C|(6*) and Cj(9) were fixed at the values 
described in the previous subsection [52Ll53| . 

For the ion-pair state, it is the quantity V^R, 9) — PJ^ 
that is asymptotically zero, where E 1 ^ = 0.0863426 P^ is 
the ion-pair threshold. This was interpolated using RP- 
RKHS interpolation with parameters m = and n = 2, 
which gives a potential with long-range form 



V\R,i 



d{9)R- 1 - C 2 (9)R- 



(18) 



The coefficient C\{9) was fixed at the Coulomb value of 
lEhdo 1 . 

The coupling potential V 1C (R, 9) has no obvious 
inverse-power form at long range. It was therefore in- 
terpolated using the ED-RKHS (exponentially decaying 
RKHS) approach 0, with n = 2 and /3 = 0.77 A. This 
gives a potential with long-range form 



V ic (R,9) = A(0)exp(-pR). 



(19) 



The value of (3 was chosen by fitting the values of the 
coupling potential at R = 10 and 12 A to decaying expo- 
nentials. 

Interpolation in 9 was carried out in a subsequent step. 
An appropriate angular form is 



V q (R,9) = J2v kq (R)<d kq (cost 



(20) 



The procedures described above produce six potential 
energy surfaces on a two-dimensional grid of geometries 
(R,9): four surfaces corresponding to covalent states 
("c"), with A' or A" symmetry and triplet or singlet 
multiplicity, one surface with ion-pair character with A' 
symmetry, and finally the non- vanishing coupling of this 
latter state with the covalent 1 A' state. We will denote 
the 6 surfaces y x K , v' iA '-, V^" , V %A ^ and V ic , 

respectively. Labels will be dropped below when not rel- 
evant to the discussion. 

For the covalent states of each spin multiplicity, inter- 
polation was carried out on sum and difference potentials, 



1 



V 2 (R,9) = ~ V A "(R,e)-V A (R,t 



V (R,9) = 1 V A "(R,t 



V A '(R,( 



(15) 
(16) 



with the difference potentials set to zero at 9 = and 
180° to suppress the slight non-degeneracy in the MRCI 



where 0fc 9 (cos( 
functions, 



Ofc g (cos( 



are normalized associated Legendre 



and q = for the sum and ion-pair potentials, 2 for the 
difference potentials and 1 for the coupling potential, q 
can thus be used as a label to distinguish the different po- 
tentials. The coefficients Vk q (R) for k = q to 9 were pro- 
jected out using Gauss-Lobatto quadrature, with weights 



Vkg(R) =X>iW>0i)e fcg (cos0 i ), (22) 



Since there are fewer coefficients than points, the result- 
ing potential function does not pass exactly through the 
potential points. However, the error for the covalent 
states in the well region is no more than 20 \iEh. 
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We thus arrive at a set of i?-dependent coefficients 
V^ 2S+1 (R), with a = c, i or ic labeling the potentials 
for the pure covalent or ion-pair states or the coupling 
between them and 25* + 1 = or 1 for singlet or triplet 
states respectively. These coefficients will be used be- 
low in evaluating the electronic potential matrix elements 
which couple collision channels in the dynamical calcula- 
tions. 

III. DYNAMICAL METHODOLOGY 
A. The basis sets 

We carry out coupled-channel calculations of the col- 
lision dynamics. The channels are labeled by quantum 
numbers that characterize the internal states of the col- 
liding partners, plus partial wave quantum numbers that 
define the way the partners approach each other. 

It is convenient to distinguish between the laboratory 
frame, whose z axis is taken to be along the direction of 
the external field (if any) and the molecule frame, whose 
z axis lies along the internuclear axis of the OH molecule 
and whose xz plane contains the triatomic system. This 
allows us to define external coordinates that fix the col- 
lision plane, and internal coordinates that describe the 
relative position of the components on it. As external co- 
ordinates, we choose the Euler angles (a, f}, 7) required to 
change from the laboratory frame to the molecule-frame; 
as internal coordinates we use the system of Jacobi coor- 
dinates (R, 8) defined above. We also define the spherical 
angles (0' ,</)') that describe the orientation of the inter- 
molecular axis in the external (laboratory) frame. 

We first focus on the covalent channels OH(7r 3 , 2 IT) + 
Rb^s 1 , 2 Si/ 2 )- The OH molecule can be described us- 
ing Hund's case (a) quantum numbers: the internal state 
is expressed in a basis set given by \sdO~)\X)\jmLo) , where 
Sd is the electronic spin of OH and a its projection on 
the internuclear axis; A is the projection of the electronic 
orbital angular momentum onto the internuclear axis; j 
is the angular momentum resulting from the electronic 
and rotational degrees of freedom, m its projection on 
the laboratory axis and lo = A + a its projection on the 
internuclear axis. The symmetric top wavefunctions that 
describe the rotation of the diatom in space are defined by 

b'ma>) = \f^D3* UJ {a,l3 ll ). At this stage A, a and lo 
are still signed quantities. However, in zero electric field, 
energy eigenstates of OH are also eigenstates of parity, 
labeled e = e or /. These labels refer to "+"or "— " sign 
taken in the combination of lo and —lo; the real parity 
is given by e(— l) J ~ Sd . To include parity, we can define 
states \s d Xu)ejm), where the bar indicates the absolute 
values of a signed quantity. Finally, we need to include 
also the nuclear spin degree of freedom: if id designates 
the nuclear spin angular momentum of the diatom, then j 
and id combine to form fd, the total angular momentum 
of the diatom, which has projection mf d on the labora- 
tory z axis. The resulting basis set that describes the 



physical states of the OH molecule is 

\s d \ue{jid)fd,nif d .) (23) 

For Rb, the electronic angular momentum is given en- 
tirely by the spin s a of the open-shell electron, which 
combines with the nuclear spin i a to form f a , the total 
angular momentum of the atom. The state of the Rb 
atom can then be expressed as 

\{s a ia)fam Sa ). (24) 

The explicit inclusion of the Sd and s a quantum numbers 
allows us to use the same notation for ion-pair channels 
OH-(tt 4 , 1 S+) + Rb + ( 1 5 ): in this case both partners 
are closed-shell, so Sd — s a — X — lo = and e = +1. 

The resulting basis set for close-coupling calculations 
on the complete system (designated Bl) has the form 

\B1) = \s d XQe(ji d )f d m fd )\(sJ a )f a m u }\LM L ), (25) 

where \LMi) denotes the partial wave degree of freedom 
and is a function of the (8',4>') coordinates considered 
above. The ket \LMl) corresponds to Ylm l (^'i 0')- Ul- 
timately, S-matrix elements for scattering are expressed 
in basis set Bl. 

The scattering Hamiltonian is block-diagonal in total 
angular momentum and total parity. The total parity is 
well defined in the basis set we have selected, and given 
by p = e(— iy~ Sd + L . It is conserved in the presence of a 
magnetic field but not an electric field. The total angular 
momentum is not conserved in the presence of either a 
magnetic or an electric field. However, the projection of 
total angular momentum on the laboratory z axis, given 
by Mj = mt d + mt + Ml, is conserved in the presence 
of an external field aligned with the laboratory z axis. 

B. Matrix elements of the potential energy 

The PES in section II are diagonal in the total elec- 
tronic spin, S, and in the states of the nuclei, rrii a and 
m,i d , with m !a and the nuclear spin projections of the 
atom and diatom in the laboratory frame. We therefore 
find it convenient to define basis sets labeled by these 
quantum numbers. This allows not only the direct calcu- 
lation of potential matrix elements, but the definition of 
some useful frame transformations & Two other basis 
sets, B2 and B3, defined with/without parity (B2p/B2w 
and B3p/B3w) are described in Appendix 2. The cor- 
responding frame transformations are defined in section 

unfa 

The calculation of the matrix elements of the potential 
energy is most direct in basis set B2w. This is based on 
Hund's case (b) quantum numbers for the molecule, and 
is given by 

\B2w) = \(s d s a )SMs)\nm n X)\i d m id )\i a m ia )\LM L ),(26) 

where n is the total angular momentum excluding 
spin, with projection A on the internuclear axis, \SMs) 
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indicates the total spin state of the electrons, and 
\ia r nia}\^d, m i d ) indicates the states of the nuclei. 

In order to relate the PES obtained in section II to 
the quantum numbers of our channels, it is convenient to 
recast the electronic wavefunctions for the covalent states 
of 2S + 1 A' and 23+1 A" in terms of functions with definite 
values of A, 

r (3) +ic) =-i(i i ^4>+*r (3) 4 / >); 

I 1(3) -lc) =±={\^K)-A^Ai>)). (27) 

We can also associate A = with the ion-pair wavefunc- 
tion (and denote it 1 1 0i) ) . Then the multipole index q 
in the potential expansion (|2UI) is viewed as an angular 



momentum transfer q = A' — A. 

The B2w basis functions do not explicitly depend on 9. 
We therefore rotate the Qkq functions |55| onto the labo- 
ratory frame, to which molecular spins and partial waves 
are ultimately referred. The functions 0fc g are propor- 
tional to renormalized spherical harmonics Ck q {d^ 0) |55| . 
for which 

C kg (e, 0) = J2 D k vq {a, /3, 7 )CW(0', tf). (28) 

V 

The potential now depends on the same angular coordi- 
nates as the B2w basis functions. Integrating and apply- 
ing the usual relationships, we obtain the matrix elements 
of the electronic potential in basis set \B2w), 



(LM L \{(s d s a )SM s \{nm n X\V\ n 'm' n X% S y' a )S'M' s )\L'M' L ) 



fok«^VH[«'][i][i'](-i) (MlK - V) i _ Ml v m , 



L k V 





k 




Hi 








n k )( ? k ) \ ¥^KV"f +1 (R), (29) 

m„ v -m' I \ X a -A' I \ 2 k ^ 



where [A] = 2 A +1, q = X' — A and v = m' n — m n . k is a 
constant whose value in the case of covalent-covalent or 
ionic-ionic matrix elements (q = 0, ±2) is 1. For covalent- 
ionic or ionic-covalent matrix elements (q = ±1), k = 
— y/l/2 or k = \/l/2 respectively. In these last two 
cases, V^ { s+1 (R) = Vlf(R). 



Our aim is to evaluate the matrix elements of the po- 
tential in basis set Bl. Basis set B3w, defined in Ap- 
pendix 2, can be considered as an intermediate step be- 
tween B2w and Bl. Starting from Ea. 1291 and changing 
basis to B3w (see Appendix 2), we obtain 



(LM i |<t 8 mU<«„mJ<»«im i J0^ 



k m 3 rn' S 



s d s a S \ ( s' d s' a S 



m s d m Sa -M s J I m 'sri m 's„ ~ M s 



x (_ 1 yn+ m '-u,-u,' + M L + m ' n -\'+j+ ] '+n+n'+s d +s d ^ ^ ^ ^ ^ ^ 



mi 1 



3 s d n \ / n k n' \ / L k V \( f s' d n' 
m -m Sd -m n J \ m n v -m' n J y -M L v M' L ) \ m' -m' Sd -m' n 

( i * ) ( - : ) ( ; s L o ) ( 5 % -y ) JW^rw- « 



The potential matrix elements in basis sets B2p and B3p 
are trivially related to the ones in B2w and B3w, respec- 
tively, requiring only the change to a parity-symmetrized 
basis set, built as a superposition of +A and —A or +u> 



and — lo vectors. Finally, the evaluation of the potential 
in basis set Bl can be easily obtained from that in B3p 
by taking the standard composition of j and s a with the 
respective nuclear angular momenta id and i a - 
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Rb(£ =2) 



OHM, fj=2)+Rb(f ,=2) 
OH(£=f,fj=l)+Rb(L=2) 



OH{e=e,fj=2)+Rb(f, =2) 
OH(E=e,fj=l)+Rb(f, =2) 



■CI 
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OH(E =f) 



Rb(f, =1) 
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OH(fd=2) 
j OH(£=e) 3l j K ; 0H (f i= i) 



OHM,fj=2)+Rb[f,=l) 
OHM,fj=l)+RblL=l) 



OH(e=e,fj=2)+Rb(f.=l) 
OH(e=e,fj=l)+Rb(f.=l) 



C2 



FIG. 8: Threshold diagram for the channels of interest, la- 
beled by quantum numbers that characterize the internal 
states of both colliding partners. These correspond to the 
ground rotational state uj = 3/2, j = 3/2 of OH. Thresholds 
corresponding to the incident channels considered in section 
II VI (CA and C2) are indicated. Hyperfine splitting in Rb and 
A-doubling and hyperfine splitting in OH are given. 



The true eigenstates of OH are linear combinations 
of functions with different values of oj, mixed by spin- 
uncoupling terms in the hamiltonian. The mixing is 
significant even for the rotational ground state: 85% of 
uj = 3/2 and 15% of uj = 1/2. We have approximated OH 
as a pure case (a) molecule for convenience in the present 
work. The fine and h ype rfine energies, taken from the 
work of Coxon et al. were associated to a unique 

set of case (a) quantum numbers \sdXuie(jid) fd.mfj). For 
Rb, the experimental values are taken from ref. |53- Fig- 
ure [S] shows the quantum numbers that characterize the 
internal states of the colliding partners for the 8 lowest 
asymptotic thresholds, corresponding to the rotational 
state lj = 3/2, j = 3/2 of OH. 

C. Solving the Schrodinger equation 

The full Hamiltonian operator can be written 



2/x 



\dR 2 R 2 



V(R) + H t m) 



where [i is the reduced mass, V(R) indicates the poten- 
tial matrix containing the electronic potential matrix el- 
ements given in the previous section and Hth is the con- 
stant thresholds matrix. By construction, any constant 
difference of energy between channels has been relegated 
to the threshold matrix so that the potential matrix ele- 
ments die off at long range. 

The coupled-channel equations that result from intro- 
ducing this Hamiltonian into the total Schrodinger equa- 
tion are propagated using the log-derivative method of 



TABLE III: Approximate number of partial waves required 
in a converged calculation for each rotational state (uJ,j)- 



j 


= 1/2 j = 3/2 j 


= 5/2 j 


= 7/2 j 


= 9/2 j 


= 11/2 


Q = 3/2 


70 


70 


60 


50 


25 


uj = 1/2 


70 70 


70 


65 


50 


40 



Johnson [53, modified to take variable step sizes keyed 
to the local de Broglie wavelength. The log-derivative 
matrix Y thus obtained is matched to spherical Bessel 
functions in the usual way to yield the scattering ma- 
trix S. Using T matrix elements (T = i(S - 1)) [60| . 
the cross section for incident energy E between internal 
states a and /3, corresponding to a beam experiment, can 
be obtained using the expression 



a a ^(E) = ^ L ""V(2^+1)(2L" + 1) 

LL'M' L L" 

X T aL0^fJL'M'( E ) T aL"Q^pL'M' L (E), (32) 



where the collision axis is chosen along the quantiza- 
tion axis (so that only Ml = contributes). The 
labels a and j3 designate sets of quantum numbers 
specifying the internal states of both colliding partners 
(s d \uje(ji d )fdmf d , (s a i Q )/ a m /a ), and k = (2^E /h 2 ) 1 / 2 is 
the wave number corresponding to incident kinetic en- 
ergy E. 

Our aim is to extract information for collisional pro- 
cesses involving OH(j = 3/2, 2 II3/2) +Rb(5s 1 , 2 S 1/2 ), in 
any of their hyperfine states, with translational energies 
in the range 10 -6 to 1 K. This range comfortably includes 
both the temperatures currently reached in buffer-gas or 
Stark deceleration experiments and the target tempera- 
tures of sympathetic cooling. A fully converged calcula- 
tion would require propagation to large values of R and 
the inclusion of a huge number of channels. Although 
both covalent and ion-pair channels should be consid- 
ered, in this pilot study we include only covalent states, 
except in Section llV Cl below. The large anisotropy of the 
surfaces included in the calculation makes it necessary to 
use a large number of OH rotational states and partial 
waves; the inclusion of rotational states up to j = 11/2 
is required (these states lie fts 550 cm -1 (oj = 3/2) and 
w 850 cm -1 (uj — 1/2) above the ground state, numbers 
comparable to the depth of the covalent potential energy 
surfaces). The number of partial waves needed for con- 
vergence for each rotational state is shown in Table ITTll 
it may be seen that fewer partial waves are required for 
higher rotational states. Unfortunately, a fully converged 
calculation was beyond our computational resources, so 
we reduced the number of channels to include only 40% 
of those in Table [TTT1 This gives cross sections accurate 
to within a factor of 2. 

We consider first the collision of atoms and molecules 
that are both in their maximally stretched states, 
(/ a ,m/ a ) = (fd,rrif d ) — (2, +2). The s-wave incident 



12 



channel has L = and Ml = 0, so corresponds to 
Mj = +4. The set of all Mj = +4 channels with a 
defined total parity p, including all allowed Ml projec- 
tions, as well as all the f a , fd, nif a and m/ d quantum 
numbers (or equivalently m, m Sa and m, n , m,i d ) contains 
23433 channels. This makes an exact calculation infcasi- 
ble. We have therefore introduced two approximations to 
reduce this number. First, the projection rrif d for chan- 
nels with j > 5/2 is fixed to its initial value; the sup- 
pressed projections increase the degeneracy and might 
split rotational Feshbach resonances, but numerical tests 
show that making this approximation does not substan- 
tially alter the overall magnitude of the cross sections 
reported here. This approximation reduces the number 
of channels to 10555. Second, for propagation at large R 
we disregard channels that are "locally closed", that is, 
whose centrifugal barrier is higher than the incident en- 
ergy in a given amount that is modified until converence. 

Even with these approximations, and the suppression 
of ion-pair channels, it is impractical to perform full cal- 
culations. However, dividing the radial solution of the 
Schrodinger equation into an inner region (R < Rq) and 
an outer region (R > Rq) makes it possible to use differ- 



ent basis sets ( "frames" ) in each of them. Frame transfor- 
mations have previously been employed for the simpler 
problem of alkaline earth+alkaline earth collisions [63] 
and for electron- molecule collisions 61] . They will be 
an essential tool for introducing hyperfine structure into 
atom-molecule and molecule-molecule collision problems 
[63| . The calculation is thus divided into two different 
steps: 



• At short range, R < Rq, the hyperfine interac- 
tion is small compared to the depths of the short- 
range potentials. We therefore represent the Hamil- 
tonian in basis set B3 (see Appendix 2), where 
the potential is diagonal in nuclear spin projec- 
tions nii a and m^. There are 8 such blocks (since 
(2id + 1)(2« + 1) = 8). We ignore elements of H t h 
that couple different pairs of rrii a and . This re- 
duces a single (8N) x (8 A) calculation to 8 (TV x N) 
calculations. At R = Ro the complete Y matrix can 
be rebuilt using the partial Y mia ,m« matrices ob- 
tained from each subset, \sdXu>ejrn)\s a m a }\LML) 



(LM L I (i a m la I (s a m Sa | {i d m id \ (s d Xujejm\Y \s d \'u'c'j'm') \i' d m' id ) \s' a m' s J \i' a m' ia ) \L'M' L ) 

= ^iai'Jmi^Jiii'Jmi^'i {LM L \{s a m s J[{sd (33) 



This Y matrix is then transformed into the asymp- 
totic basis set Bl. We have found that this frame 
transform provides a very accurate way to include 
the Rb-OH hyperfine structure in reduced calcula- 
tions using only covalent channels. Moreover, ow- 
ing to the depth of the short-range potentials, Y 
is weakly dependent on energy and can be interpo- 
lated in the inner region. 

• At long range, R > Rq, the Y matrix, expressed 
already in Bl, is propagated to large distances 
to obtain the S matrix. We invoke an alterna- 
tive approximation: The coupling between different 
asymptotic rotational and fine-structure states di- 
minishes at longer distances and can be neglected. 
Thus the subset corresponding to the ground rota- 
tional diatomic state (j = 3/2, uj = 3/2) can be 
propagated by itself to asymptotic distances. We 
reinstate the full hyperfine Hamiltonian in this re- 
gion. 

It is worth noting that basis set B2, defined in Ap- 
pendix 2, could also be the basis for another frame trans- 
formation. Approximate decoupling of singlet and triplet 
channels in the inner region would allow a partition of the 
numerical effort into two smaller groups of channels, and 



the introduction of the ionic channels in the singlet one. 

IV. SCATTERING CROSS SECTIONS 

We have calculated elastic and state-resolved inelastic 
cross sections for two different incident channels for col- 
lisions of Rb atoms with OH molecules. These are shown 
as CI and C2 in Figure [S] Although we do not include 
the effects of external fields explicitly in this work, we 
consider states in which both partners can be in weak- 
field-seeking states in a magnetic field, and thus magnet- 
ically trappable. The OH hyperfine states that can be 
trapped at laboratory magnetic fields are (fd — 2, mf d = 
+2, +1, 0) and (fd = l,rrif d = +1), while the correspond- 
ing states for Rb are (f a — 2,m/ a = +2, +1,0); the Rb 
state (f a — \,m$ a — —1) is trappable for fields smaller 
that » 1250 gauss. 

In the first case, designated CI in Figure|Sl both part- 
ners are in maximally stretched states: OH(e = f, fd — 
2,mf d = 2) + Rb(/ Q = 2,m/ Q = 2). This case corre- 
sponds to the highest threshold correlating with OH in 
its ground rotational state. In this case, OH is also elec- 
trostatically trappable. The second case, designated C2 
in FigureEl correlates with the lowest asymptotic thresh- 
old in the absence of an external field: OH(e = e, fg = 
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R(a.u.) 

FIG. 9: Adiabatic curves correlating with the lower rotational 
states for the collision in maximally stretched states of both 
partners (p = +l,Mj = +4): for R < 17 ao hyperfine inter- 
action is partially neglected; for long range R > 17 ao, only 
channels with OH in its ground rotational state are consid- 
ered and full hyperfine is included (inset gives a more detailed 
image of the latter regime). 

l,m/ d = 1) + Rb(/ a = l,m/ a = —1). Both partners are 
again magnetically trappable, although in a field they 
will no longer be the lowest energy states. 

Figure shows selected adiabatic curves correlating 
with the lower rotational states for the collision with both 
partners in maximally stretched states (p = +l,Mj = 
+4). In our calculations we take Rq = 17 a . As de- 
scribed above, for R < Rq the hyperfine interaction is 
partially neglected. For R> R , only hyperfine channels 
with OH in its ground rotational state are included. 

A. General behavior: total cross sections 



expected in this region. Rather, the values of threshold 
cross sections are strongly subject to details of the PES, 
and are typically only uncovered by experiments. 

At higher energies, above 1CP 2 K, where many par- 
tial waves contribute (I > 4 for CI and I > 5 for C2), 
the behavior of the cross sections changes and inelastic 
processes are well described by a semiclassical Langevin 
capture model 33, 36], 

CTLangevin (E) = 3?T ( — J . (34) 

This cross section is also plotted in Figure ITU1 (points). 
The Langevin expression arises as a limit of the exact 
quantum expression in Eq. 1321 with the usual assump- 
tions: the impact parameter takes values in a continu- 
ous range, and the height of the centrifugal barrier, de- 
termined using only long-range behavior, determines the 
number of partial waves that contribute for a given en- 
ergy. Similar behavior has been observed previously in 
cold collisions, such as M + Mi with M an alkali metal. 
For Li + Li2, I > 3 was found to be a sufficient for the 
cross sections to exhibit Langevin behavior |30|. 

As can be seen in Figure 111)1 the Langevin limit re- 
produces the general trend across the entire semiclassi- 
cal energy range. In a Hund's case (a) system like OH, 
where the electron spin is strongly tied to the intermolec- 
ular axis, the highly anisotropic PES might be expected 
to disrupt the spin orientation relative to the laboratory 
frame completely. As a consequence, inelastic processes 
are expected to be very likely and the Langevin model 
should describe well the behavior of Rb-OH and similar 
systems. A similar upper limit for the elastic cross sec- 
tion is given by four times the inelastic one. It is easy to 
verify that, if the inelastic cross section reaches its maxi- 
mum value, the elastic and inelastic contributions to the 
cross section must be equal. This behavior is also seen 
in Figure ITU1 

The cross sections for incident channel C2 are quite 
different from those for incident channel CI. For C2, the 
cross sections are highly structured, exhibiting a large 
number of Feshbach resonances. Since the atom and the 
diatom are both in their lowest-energy state, there are 
plenty of higher-lying hyperfine states to resonate with. 
However, the pronounced minimum in the elastic cross 
section at E ~ 10~ 4 K is the consequence of a near-zero 
s-wave phase shift. 



We begin by showing in Figure EH the total cross sec- 
tions for incident channels CI and C2. A brief description 
of these has been reported previously [3i| ■ Below 10~ 4 K 
for incident channel CI, or 10 -5 K for C2, the Wigner 
threshold law applies. Namely, as the energy goes to 
zero, cross sections corresponding to elastic and isoen- 
ergetic processes approach a constant value, while those 
for exoergic processes vary as rapidly exceeding 

elastic cross sections. No quantitative predictive power is 



B. Detailed picture: partial cross sections 

Some additional insight into the collision process can 
be obtained by examining the partial cross sections to 
various final states. However, there are many of these. 
Starting in incident channel CI, there are 128 possible 
outcomes, counting the hyperfine states of both Rb and 
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FIG. 10: Elastic and total inelastic cross sections for the cases 
analyzed: Panel (a) corresponds to OH(e = f,fd — 2, m/ d 
•) + Rb(/ a = 2,m fa = 2) incident channel (CI). Panel 
(b) corresponds to OH(e = e, f d — l,m/ d = 1) + Rb(/ a = 
b m /a — — 1) incident channel (C2). The points indicate the 
Langevin cross section. 



OH, plus the lambda doublet of OH. To simplify this 
information, we first break the cross sections into four 
classes: elastic scattering, scattering in which only the 
Rb state changes, scattering in which only the OH state 
changes, and scattering in which both change. 

These four possibilities are shown in Figure lTTl for both 
incident channels, CI and C2. In general, all four pro- 
cesses are likely to occur. This further attests to the 
complete disruption of the spin during a collision. Nev- 
ertheless, at the highest energies probed (where results 
are less sensitive to potential details) there is a definite 
propensity for the OH molecule to change its internal 
state more readily than the Rb atom, at least when only 
one of them changes. This is probably a consequence of 
the spherical symmetry of the Rb atom, whereby its elec- 
tronic spin is indifferent to its orientation. By contrast, 
the electronic angular momentum of OH is strongly cou- 
pled to the molecular axis, and will follow its changes in 
orientation due to the anisotropies in the interaction. 

A more detailed understanding can be obtained by 



FIG. 11: Total cross sections for 4 types of process are con- 
sidered: those where only OH partner changes internal state, 
those where only Rb changes, those where both partners suf- 
fer inelastic change and those where neither change. Panel 
(a) corresponds to OH(e = /, fd = 2, mf d =2) + Rb(/ a = 
2, m/ = 2) incident channel (CI). Panel (b) corresponds to 
OH(e = e, fd — l,mf d — 1) + Rb(/ a = l,m,f a = —1) incident 
channel (C2). 



considering Eq. 1301 For Rb the hyperfine projection is 
given by TO/ a = m Sa + rrii a , whereas for OH it is given 
by = m n + m Sd + rrn d . The nuclear spin projec- 
tions rrii a and are untouched by potential energy 
couplings. The potential conserves Ms = m Sa + m Sd , so 
that if the Rb electronic spin m Sa changes, so will the 
OH electronic spin m Sd . On the other hand, OH can also 
change the projection m n of the rotational angular mo- 
mentum n, which is absent in Rb. Thus OH has more 
opportunities to change its internal state than does Rb, 
and this is reflected in the propensities in Figure ITT1 

We first consider incident channel CI, with both colli- 
sion partners initially in their maximally stretched states. 
Some of the main outcomes are shown in Fig^J In the 
absence of an anisotropic interaction, and remembering 
that Mj = rrif a +rrif d + Ml is conserved, the sum of the 
spin projections Mp = mt + rrif d would be conserved. 
Thus inelastic collisions would be impossible. This is be- 
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FIG. 12: Dominant partial cross sections for OH(e = /, fd 
2,m fd =2) + Rb(/ a = 2,m fa = 2) incident channel (CI). 



cause the projection of m / d could not be lowered without 
raising mf a , but mj a could not be raised further. How- 
ever, the anisotropic PES in Rb-OH allows such changes 
quite readily. 

There remains, however, a propensity for collisions 
with small values of Am/ to be more likely, as seen in 
Figure ED Part (a) shows cross sections that change 
m f d for OH without changing m t a for Rb, while part (b) 
shows those that change m f a without changing m t d . As 
noted above, Rb appears to be more reluctant to change 
its projection. In fact, since the potential is diagonal 
in the states of the nuclei, only consecutive values of 
mf a are coupled in first order. On the other hand, first- 
order coupling exists between many different values of 
m f d . Finally, a decrease in the probability of processes 
with increasing AMp can be related to the diminution 
of anisotropic terms in the potential when increasing the 
angular momentum transfer (k in Eq. 130(1 . 

We now consider incident channel C2, OH(e = e, fd = 
l,m/„ = 1) + Rb(/ a = l,mf a = -1). This is the 
lowest threshold in the absence of a field. Only three 
channels, degenerate with the initial channel, arc possi- 
ble outcomes at very low energies: |e = e, 1, +1) |1, — 1) 
\e = e,l,0)|l,0), and \e = e, 1, -1)|1, +1). Since these 
states have the same value of Mp = mf d +mf a as the ini- 
tial channel, the processes can occur by ordinary spin ex- 
change with no centrifugal barrier. Partial cross sections 
for these three channels are shown in Figure 1141 (a) over 
the entire energy range. Intriguingly, inelastic (state- 
changing) collisions seem to be somewhat suppressed 
relative to elastic scattering. Suppressed spin-exchange 
rates would presumably require delicate cancellation be- 
tween singlet and triplet phase shifts |7^|. That such a 
cancellation occurs in a highly multichannel process is 
somewhat unexpected. However, examining the matrix 
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FIG. 13: Panel (a) ((b)): inelastic cross sections for mf d 
( m f a ) changing collisions. The higher the change, the lower 
the cross section. They correspond to OH(e = /, /<j = 
2,m /d = 2) + Rb(/ a = 2,m u = 2) incident channel (CI). 



elements of the potential reveals that there is no direct 
coupling between the initial state |e = e, 1, +1) |1, — 1) 
and final state |e = e, 1, — 1)|1, +1). Transitions via po- 
tential couplings are therefore a second-order process re- 
quiring the mediation of other channels. 

Many other exit channels are possible, once energy and 
angular momentum considerations permit them. Cross 
sections for several such processes are shown in Figure 
1141 (b). For example, the channels OH(e = e,l,0) + 
Rb(l,+1) and OH(e = e,l,0) + Rb(l,-1) are not con- 
nected to the initial channel by spin exchange, since Mp 
changes by —1. In this case, angular momentum shunts 
from the molecule into the partial-wave degree of free- 
dom, necessitating an L = 1 partial wave in the exit 
channel. Therefore, this process is suppressed for ener- 
gies below the p-wave centrifugal barrier, whose height is 
1.6 mK. 

We conclude this subsection by stressing the vital im- 
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threshold. Only after both channels are open, at higher 
energies, do the cross sections become almost identical. 
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FIG. 14: Some partial cross sections for C2 incident chan- 
nel: OH(e = ej d = l,m fd = 1) + Rb(/„ = l,m /a = -1). 
On panel (a) elastic cross section together with partial cross 
sections for two degenerate processes with the initial chan- 
nel, opened at zero collision energy, are shown. The processes 
shown on panel (b) are also degenerate but closed by the cen- 
trifugal barrier at zero collision energy. 



portance of including the hyperfine structure in these 
calculations. Figure El shows partial cross sections 
for incident channel CI scattering into pairs of chan- 
nels which differ only in the small hyperfine splitting 
of OH (w 3 mK). In one example (black and red solid 
curves), the parity of OH changes from / to e, leading 
to the final channels OH(e = e,fd = l,2,m/ d = 0) + 
Rb(/ a = 2, m/ a =0). These channels, distinguished only 
by their hyperfine quantum number fd, are almost iden- 
tical at high energies, but quite different at low ener- 
gies. As a second example, consider the final channels 
OH(e = /,/„= l,2,m fd = 0) + Rb(/„ - 2,m fa = +2), 
(green and orange dashed) , which preserve the initial par- 
ity. The process with fd = 1 is exothermic, while the one 
with fd = 2 requires the opening of the partial wave 
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FIG. 15: Partial cross sections for two pairs of processes corre- 
sponding to incident channel CI (OH(e = /, fd — 2, rrif d = 2) 
+ Rb(/ a = 2,m/ a =2)) and final channels differing only in 
the fd quantum number. Their values converge at high kinetic 
energies but are very different in the cold regime. 



C. The harpooning process 

We have not yet fully incorporated the ion-pair chan- 
nel in our calculations, but it is instructive to assess its 
influence. In the harpooning model it is usually under- 
stood that, if the electron transfer takes place at long 
range, large collision cross sections result j64(. However, 
for Rb-OH the crossing point i?o ~ 6 A gives a geometric 
cross section of only a = 4itRq ps 4.5 • 1CP 14 cm 2 . This is 
substantially smaller than the quantum-mechanical cross 
sections we have estimated. Thus while harpooning may 
enhance cross sections substantially at thermal energies, 
it is likely to have less overall influence in cold collisions. 

The inelastic results above place cross sections near 
the semiclassical Langevin upper limit. The harpooning 
mechanism is unlikely to increase their magnitude, but 
could change their details. For example, it might be ex- 
pected that the harpooning mechanism would distribute 
the total probability more evenly among the different spin 
orientations. 

We have explored the influence of the harpooning 
mechanism in a reduced calculation where the hyper- 
fine structure is neglected. This makes the calcula- 
tion tractable even with the ion-pair channel included. 
The resulting elastic and total inelastic cross sections, 
for incident channel OH(e = f,j = 3/2, m — 3/2) + 
Rb(s a = 1/2, m Sa — 1/2) (maximally stretched in the 
absence of nuclear degrees of freedom) are shown in Fig- 
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urc ll6l They are compared to the results obtained with- 
out the ion-pair channels. In the semiclassical region, 
E > 1 mK, the general order of magnitude of the cross 
sections is preserved, although the detailed features are 
not. In the Wigner regime, by contrast, including the 
ion-pair channel makes a quite significant change, be- 
cause completely different phase shifts are generated in 
low-lying partial waves. We stress again, however, that 
we do not expect this or any model to have predictive 
power for the values of low-energy cross sections until the 
ambiguities in absolute scattering lengths are resolved by 
experiments. 
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FIG. 16: Effect of the inclusion of ion-pair channels in a re- 
duced (neglecting hyperfine structure) calculation. Langevin 
cross section is also shown for comparison. 



V. CONCLUSIONS AND PROSPECTS 

Based on the results above, we can draw several general 
conclusions concerning both the feasibility of observing 
Rb-OH collisions experimentally in a beam experiment 
and the possibility of sympathetic cooling of OH using 
ultracold Rb. 

We can assert with some confidence that Rb-OH cross 
sections at energies of tens of mK, typical of Stark de- 
celerators, will be on the order of 10~ 13 cm 2 . This is 
probably large enough to make the collisions observable. 
In round numbers, consider a Rb MOT with density 
^Rb = 10 10 cm -3 , and an OH packet emerging from a 
Stark decelerator at 10 m/s, with density noH = 10 7 
cm~ 3 . Further assume that the Rb is stored in an elon- 
gated MOT that allows an interaction distance of 1 cm 
as the OH passes through. In this case, a cross section 
a = 10 -13 cm 2 implies that approximately 10~ 3 of the 
OH molecules are scattered out of the packet. This quan- 
tity, while small, should be observable in repeated shots 



of the experiment 65] . 

As for sympathetic cooling OH using cold Rb, it seems 
unlikely that this will be possible for species in their 
stretched states as in channel CI. The inelastic rates are 
almost always comparable to, or greater than, the elastic 
rates. The situation becomes even worse as the tem- 
perature drops and cross sections for exoergic processes 
diverge. Every collision event might serve to thermalize 
the OH gas, but is equally likely to remove the molecule 
from the trap altogether, or contribute to heating the 
gas. 

For collision partners in the low-energy states of inci- 
dent channel C2, the situation is not as bleak, at least 
in the low-energy limit, since inelastic cross sections do 
not diverge and may be fairly small. However, we have 
focused on the weak-magnetic-field seeking state of OH, 
which will not be the lowest-energy state in a magnetic 
field. Thus exoergic processes again appear, and inelas- 
tic rates will again become unacceptably large. Under 
these circumstances, rather than "sympathetic cooling," 
the gas would exhibit "simply pathetic cooling." 

From these considerations, it seems that the only way 
to guarantee inelastic rates sufficiently low to afford sym- 
pathetic cooling would be to remove all exoergic inelas- 
tic channels altogether. Thus both species, atom and 
molecule, should be trapped in their absolute grou nd 
states, using optical or microwave dipole traps |6|| or 
possibly an alternating current trap [67j. Sympathetic 
cooling would then be forced due to the absence of any 
possible outcome other than elastic and endoergic pro- 
cesses. The latter might produce a certain population of 
hyperfine-excited OH molecules, which could then give 
up their internal energy again and contribute towards 
heating. However, the molecules are far more likely to 
collide with atoms than with other molecules, so that 
this energy would eventually be carried away as the Rb 
is cooled. 

Further aspects of the Rb-OH collision problem need 
further work. Foremost among these is the possibility 
that inelastic rates could be reduced or controlled by 
applying electric and/or magnetic fields. This may re- 
sult partly from the simple act of moving Feshbach res- 
onances to different energies, or from altering the effec- 
tive coupling between incident and final channels, as has 
been hinted at previously 13, 68] . In addition, the influ- 
ence of the ion-pair channel may be significant. There is 
the possibility, for example, that the highly polar RbOH 
molecule might be produced by absorption followed by 
either spontaneous or stimulated emission |37| . Finally, 
new phenomena are possible, such as analogues of the 
field-linked states that have been predicted in dipolc- 
dipole collisions [69l l70|] . 
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Appendix 1 

As described above, we encountered difficulties with 
non-degeneracies between the 1 A' and 1 A" components of 
the 1 II state at linear geometries. These could have been 
avoided by carrying out the linear calculations in 
rather than C s symmetry, but C s symmetry was essential 
to avoid discontinuities at 8 = and 180°. 

To understand the non-degeneracies it is essential to 
understand the procedure used to generate the basis set 
used in the CISD calculation. 

1. The molecular orbital basis set is partitioned into 
internal and external sets, and the internal orbitals 
are in turn partitioned into closed and active sets. 

2. A reference space is generated, including all pos- 
sible arrangements of the N available electrons 
among the active orbitals that give states of the 
specified symmetry and spin multiplicity. The 
closed orbitals are doubly occupied in all reference 
configurations. 

3. A set of (7V-2)-electron states is generated by all 
possible 2-electron annihilations from the reference 
states (with no symmetry constraints). Some of 
the closed orbitals may be designated as core or- 
bitals, in which case annihilations from them arc 
not included. 

4. A set of (TV-l)-electron states is generated by all 
possible 1-electron additions to internal orbitals of 
the (7V-2)-electron states (with no symmetry con- 
straints). 

5. The final CI basis includes states of 3 different 
classes: 

• Internal states are generated by all possi- 
ble 1-electron additions to the (TV-l)-electron 
states, with the extra electron in an internal 
orbital, that give states of the specified sym- 
metry and spin multiplicity; 

• Singly external states are generated by all 
possible 1-electron additions to the (N-l)- 
electron states, with the extra electron in an 
external orbital, that give states of the speci- 
fied symmetry and spin multiplicity; 

• Internally contracted doubly external states 
are generated by 2-electron excitations into 
the external space from a reference function 
obtained by solving a small CI problem in the 
reference space. 

The non-degeneracies between the 1 A' and 1 A" com- 
ponents of the 1 II state arise from two different sources. 



First, there is a part of the non-degeneracy due to the 
internal and singly external configurations. At linear ge- 
ometries there are both S-type and II-type reference con- 
figurations. Both types are included for A' symmetry, 
but only the II-type reference configurations are included 
for A" symmetry. In the A' calculation, therefore, there 
are additional II-type internal and singly external config- 
urations that arise from S-type reference configurations. 
This effect could in principle be suppressed in MOLPRO 
by including reference configurations of both symmetries 
in both calculations at all geometries, but this was pro- 
hibitively expensive in computer time. In addition, as 
described below, it is responsible for only about 20% of 
the total non-degeneracy. 

Secondly, there is a part of the non-degeneracy due to 
the (contracted) doubly external configurations. It must 
be remembered that we wish to transform the I A 1 and 
2A' singlet states into a diabatic representation. To do 
this meaningfully, we need them to be calculated using 
identical basis sets. If we do the MRCI calculations sep- 
arately, this is not true: even if the reference space is 
the same, each state produces a different set of internally 
contracted doubly external states. It is thus necessary to 
calculate the two X A' states in the same MRCI block, so 
that the basis set contains both sets of contracted func- 
tions. However, this results in a larger and more flexible 
basis set of contracted doubly external functions than is 
generated in the 1 A" case. 

It is helpful to document the magnitude of the non- 
degeneracies in various possible calculations. For RbOH 
at 9 = and R = 12 A: 

• If calculations are carried out in Ci v symmetry, the 
ion-pair and covalent states have different symme- 
try and the non-degeneracy is not present; 

• For C s symmetry with the two 1 A' states calcu- 
lated in a single MRCI block, both sources of non- 
degeneracy are present and the non-degeneracy is 
22.78 fiE h ; 

• For C s symmetry with the two 1 A' states calculated 
in separate MRCI blocks, only the first source of 
non-degeneracy is present and the non-degeneracy 
is 4.08^; 

• For C s symmetry with the two 1 A' states calcu- 
lated in a single MRCI block, but reference con- 
figurations of both symmetries included, only the 
second source of non-degeneracy is present and the 
non-degeneracy is 18.70 [lEh; 

• For C s symmetry with the two 1 A' states calcu- 
lated in separate MRCI blocks, with reference con- 
figurations of both symmetries included, there is no 
non-degeneracy. 
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Appendix 2 

The interaction potential does not involve the spins of 
the nuclei and is diagonal in the total electronic spin S. 
This makes it expedient to define two additional basis 
sets. 

Basis set B2 is based on Hund's case (b) quantum num- 
bers for the molecule: the orbital angular momentum of 
the electron and the rotational angular momentum cou- 
ple to form n, with A projection on the internuclear axis. 
n coupled to the spin of the electron Sd would give us j 
for the diatomic fragment, as in basis set Bl. Instead, in 
B2 we couple s d with s a to obtain the total spin of the 
electrons, described by kets \SMs). However, we leave 
aside the states of the nuclei, which are represented by 
\ia m i a )\id m i d ) j with m ia and m id being the projections 
of the nuclear spin of the atom and diatom in the labo- 
ratory frame. Basis set B2w is then given by 

\B2w) = \{s d s a )SMs)\nm n \)\i d m id )\iami a )\LM L ),{2,h) 

where \nm n \) represents ^ ^ ^^(a, P, 7)- The par- 
ity operator E* acts on these states as follows: 

E*\{sdS a )SMs)\nm n \)\i d mi d )\i a m ia )\LM L ) = 
{-l) n+L \{s d s a )SM s )\nm n - \)\i d m id )\i a m ia )\LM L ), 

(36) 



which allows us to construct combinations of define par- 
ity, 

\B2p) = \(sdS a )SM s )\nm n \£)\i d m td )\i a m la )\LM L ).(Z7) 

The third basis set, B3, is intermediate between Bl 
and B2. Here the molecule is again labeled by quantum 
numbers corresponding to a Hund's case (a) molecule, as 
in Bl, 

\B3w) = \sd<7)\\)\jmLu)\idm id )\s a m Sa )\i a mi a )\LM L )(38) 

However, the part describing the OH fragment is taken 
as \jmui)\\)\sd<j) , that is, with signed values of A, a, and 
u>. In this case, the parity operator E* acts as follows, 

E* \sdcr) I A) \jmu) \idm id )\s a m Sa )\i a m ia )\LM L ) 

x \s d -a)\ - X)\jm - w)\i d m id )\s a m Sa )\i a mi a )\LM L ), 

(39) 

so that we can again build combinations of defined parity, 

\B3p) = \s d \Qejm)\i d m ld )\s a m Sa )\i a m ia )\LM L ). (40) 

The change from B2w to B3w, leaving aside the partial 
wave and nuclear spin states, is given by 



\s d a) \\}\jmuj)\s a m Sa ) 

yyyy[ s d s a S \ / s d n j \ / s d n j \ 
2-^2-^2-^2^, \ m Sd m Sa -M s ) \ m Sd m n -m J \ -a -X uj I 

m 3d n S Ms v 

x (_i)(m—+--.+^)| nmn A)|(a (1 * )5Ms>. 
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